Order-disorder transition in a model with two symmetric absorbing states 
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We study a model of two-dimensional interacting monomers which has two symmetric absorbing 
states and exhibits two kinds of phase transition; one is an order-disorder transition and the other is 
an absorbing phase transition. Our focus is around the order-disorder transition, and we investigate 
whether this transition is described by the critical exponents of the two-dimensional Ising model. 
By analyzing the relaxation dynamics of "staggered magnetization," the finite-size scaling, and the 
behavior of the magnetization in the presence of a symmetry-breaking field, we show that this model 
should belong to the Ising universality class. Our results along with the universality hypothesis 
support the idea that the order-disorder transition in two-dimensional models with two symmetric 
absorbing states is of the Ising universality class, contrary to the recent claim [K. Nam et al., J. 
Stat. Mech.: Theory Exp. (2011) L06001]. Furthermore, we illustrate that the Binder cumulant 
could be a misleading guide to the critical point in these systems. 

PACS numbers: 05.70.Ln, 64.60.Ht, 64.60.Cn 
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I. INTRODUCTION 

During recent decades, intense theoretical efforts have 
been devoted to classifying absorbing phase transitions 
(APTs) and, as a result, several universality classes have 
been found (for a review, see, e.g., Refs. pHf). Although 
a firmly established classification principle is still desired, 
symmetry is unequivocally expected to play an impor- 
tant role in determining universality classes. On the 
one hand, if a model which has a single absorbing state 
and which does not have extra symmetry or conserva- 
tion laws exhibits an APT, it is argued that this model 
should belong to the directed percolation (DP) univer- 
sality class On the other, many systems with two 
symmetric (sets of) absorbing states are known to form 
another universality class. Such examples with Zi sym- 
metry are the probabilistic cellular automata model Q, 
the interacting monomer-dimer (IMD) model @, the 
nonequilibrium kinetic Ising model Q, and the interact- 
ing monomer-monomer model with infinitely many ab- 
sorbing states [lOj, to cite only a few. 

Although symmetry seems important, there are also 
some systems with Z 2 symmetry which do not share criti- 
cality with the above-mentioned models. Such exceptions 
can be found in Ref. [ll[. Hence, symmetry alone is not 
sufficient to determine the universality class and further 
studies are necessary to determine a guiding principle in 
terms of symmetry. 

The starting point to develop a principle would be 
construction of a (coarse-grained) field theory for each 
universality class, followed by a renormalization group 
analysis. The connection between DP and the reggeon 
field theory has been clarified long time ago [12[ (see 
also [a, la| ) . Endeavors have also been made to formulate 
a field theory for systems with ^-symmetric absorbing 
states. Cardy and Tauber [l3|, [3 developed the field 
theory for the branching annihilating random walks with 
even numbers of offspring (B AWE) [15| , which has mod- 
2 conservation of particle number. Since Zi symmetry 



entails the mod-2 conservation of domain walls in one 
dimension, the BAWE in one dimension belongs to the 
same universality class as systems with two symmetric 
absorbing states. The mod-2 conservation, however, has 
nothing to do with the Zi symmetry in higher dimen- 
sions. In particular, the BAWE in higher dimensions 
shows trivial transitions [lU Il4| . which is not the case 
for two-dimensional models with Zi symmetry [16Hl9l |. 
In this regard, the field theory of the BAWE cannot be a 
coarse-grained description for models with Zi symmetry 
in d dimensions. Thus, there was a theoretical request to 
develop a field theory with Z2 symmetry in higher dimen- 
sions and, as a response, a phenomenological Langevin 
equation was introduced [13] (see also Ref. [13] for an 
analysis of the corresponding field theory by the nonper- 
turbative renormalization group method) . 

Recently, however, this phenomenological Langevin 
equation description has been challenged The 
Langevin equation predicts that there are in general two 
transitions in two dimensions; an order-disorder transi- 
tion which is concomitant with the Zi symmetry break- 
ing (SB) and an APT (see Ref. pi] for the first obser- 
vation of two transitions in a microscopic model). Nu- 
merical analysis of the Langevin equation revealed [l7j 
that the Zi SB transition is of the Ising class and the 
APT is of the DP class. Although a Monte Carlo sim- 
ulation of the two-dimensional IMD model also found a 
Zi SB transition followed by an APT, it was claimed 
that the critical behavior for the Z2 SB transition is not 
of the Ising class [l9j]. In fact, no Monte Carlo simula- 
tions studies up to now, to our knowledge, have clearly 
shown that the Zi SB occurring in a system with two 
symmetric absorbing states is described by the Ising crit- 
ical exponents, which was the motivation of Ref. [l9j. If 
the claim in Ref. turns out to be true, a different 
coarse-grained description from that suggested is called 
for. Even more seriously, the conclusion in Ref. [13 ] ques- 
tions the validity of the theory that any continuous Zi SB 
transition between ordered and disordered phases should 
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be described by the scalar (j) 1 theory [21|, which was the 
motivation to introduce the model- A type (according to 
the Hohenberg-Halperin classification scheme [13]) inter- 
action to the Langevin equation [l7j . 

Hence, it is necessary to study the Z2 SB transition 
exhibited by a two-dimensional model with two sym- 
metric absorbing states more extensively to make a firm 
conclusion concerning the universality class. In this pa- 
per, we thoroughly investigate the SB transition. The 
model studied in this paper will be called the two dimen- 
sional interacting monomers (2DIM) model which is a 
two-dimensional version of the model studied in Ref. [23[ . 

The paper is organized as follows: After introducing 
a model and appropriate order parameters in Sec. [Til wc 
present numerical analysis of the SB transition in Sec. Mil 
Section IIVI discusses the claim in Ref. [19( , studies the 
absorbing phase transition to confirm the DP transition, 
and then summarizes the paper. 



II. MODEL 

The 2DIM model is defined on a square lattice with 
size L 2 [L is assumed to be even). Every site is indexed 
by a two dimensional vector x = with integer com- 
ponents i and j = 0, 1, . . . , L— 1). Periodic boundary 
conditions are assumed. For later purposes, the lattice is 
subdivided into two sublattices E and O. The sublattice 
E (O) is defined as a set of sites x = with i + j 

even (odd). Each site is either occupied by a monomer 
or vacant. Two monomers are not allowed to occupy a 
single site at the same time. Each site is given a state 
variable a x which takes the value 1 (0) if site x is occu- 
pied (vacant). A configuration is characterized by state 
variables at all sites. 

The dynamic rules are as follows: A monomer attempts 
to adsorb on a randomly chosen vacant site (called a tar- 
get site). Depending on the number n of occupied nearest 
neighbors of the target site, the fate of the monomer will 
be different. If n — 0, the monomer adsorbs with rate 
1. On the other hand, if n ^ 0, the monomer adsorbs 
with rate nA n , but the adsorbed monomer immediately 
forms a dimer with a randomly chosen monomer among n 
monomers on the nearest neighbor sites, and the dimer is 
desorbed in no time. Effectively, an adsorption event on a 
vacant site with occupied nearest neighbors removes one 
monomer from the lattice. If all nearest neighbors of the 
target site are occupied, that is, if n — 4, a monomer is 
not allowed to adsorb on the target site, which amounts 
to setting A4 to zero. Since A4 = 0, any configuration 
with all vacant sites surrounded by monomers is an ab- 
sorbing state. 

To write the master equation for the above dynamics 
in a succinct way, we introduce a mathematical notation; 
for any configuration C with the state variable a y for 
every site y, C x stands for the configuration obtained by 
changing the state variable at site x to 1 — a x and by 
keeping all other state variables the same as in C. Using 



this notation, the master equation can be written as 
j t P(C;t) = ]T [W cc ^P(C x ;t) - W c „cP{C;t)\ , (1) 

X 

with the transition rate 
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W Ca! c = ^,0^,0 + Sa*,i ' ^kSa y ,o5 nyik , (2) 

v fe=i 

where ' means the sum over nearest neighbor vectors 
y of site x, a x (a y ) is the state variable at site x (y) 
in the configuration C, n x (n y ) means the number of 
occupied nearest neighbors of site x (y) in C, and 5 is the 
Kronecker delta symbol. By observing that (C x ) x = C. 
one can easily find the transition rate Wcc^ ■ 

To simulate the master equation, we have used the fol- 
lowing algorithm. First, we make a list of vacant sites 
with at least one vacant nearest neighbor. For conve- 
nience, we will refer to such a vacant site as an active 
site. Assume that there are N t active sites at time t. 
A target site out of Nt active sites is selected at random 
with equal probability. If all nearest neighbors of the tar- 
get site are vacant, it becomes occupied with probability 
At which is defined as 



max(l, Ai, 2A2, 3A3) 

If the target site has n occupied nearest neighbors (n = 1, 
2, or 3), a configuration change can occur with probabil- 
ity n\ n At. If a change is destined, one monomer out 
of n is selected with equal probability and it is removed 
from the system, which mimics the dimer desorption ex- 
plained above. After the above attempt, time increases 
by At/Nt and the list of active sites is updated in an 
appropriate way. We repeat the above procedure until t 
exceeds a preassigned maximum observation time or no 
active site exists in the system. For convenience, we set 
Ai = 2A 2 = 3A 3 = A in what follows and study phase 
transitions by tuning A. 

Now we will specify the initial condition studied in 
this paper. At t = 0, the sublattice O is empty, but a 
site in the sublattice E is occupied with probability mo 
(0 < too < 1). With this initial condition, there are only 
two absorbing states; the sublattice E is fully occupied 
and the sublattice O is empty, or vice versa. In this sense, 
absorbing states have perfect 'anti-ferromagnetic' order. 
The absorbing state with the sublattice E (O) filled with 
monomers will be called the even (odd) absorbing state. 

Since we expect two transitions (an APT and a Z2 SB 
transition), two quantities that are respectively called the 
density of active sites and the 'staggered magnetization' 
will be measured during simulations. 

The density of active sites at time t is defined as 

<f>(t,L) = ^$^„.,o(l -««.,*), (4) 
x 

where n x is the number of occupied nearest neighbors of 
site a; in a configuration at time t. If the system is in one 
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FIG. 1. Phase diagram of the 2DIM model. A c is the tran- 
sition point of the order-disorder transition and X a is that of 
the absorbing phase transition. 



of the two absorbing states at time t, cf>(t, L) is obviously 
zero. We define the (averaged) density of active sites in 
the thermodynamic limit as 



Pit) 



lim <0(i,L)>, 

L—too 



(5) 



where (. . .) stands for the average over all independent 
realizations. 

The staggered magnetization (SM) is defined as 



MU.D = -L 



E-E 

.xeE xeO. 



[2a x (t) - 1] 



(6) 



If the system is in the even (odd) absorbing state, the 
SM is M{t,L) = 1 (-1). The (averaged) SM at time t 
in the thermodynamic limit is defined as 



m(f) 



lim (M(t,L)). 

L— >oo 



Since the initial condition gives 
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(7) 

(8) 
(9) 



we get 

(M(0,L)) = m o , (0(0, L)) 



1 



2-o (1 



il) . (10) 



Note that with the above initial condition, m(i), which is 
defined in the thermodynamic limit, will remain positive 
for all finite t if < tuq < 1 and the order-disorder phase 
transition is defined by the infinite-time limit of m(i) 
such that 

(nonzero, ordered phase, 
VJl = hm m(t) — < (11) 
*->°o I 0, disordered phase. 

As we will see later, the 2DIM model exhibits two 
transitions: an order-disorder transition occurring at 
A = A c and an absorbing phase transition occurring at 
A = A Q < A c . The schematic phase diagram is depicted 
in Fig. HJ 
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FIG. 2. (Color online) Plots of m(t)t pnvz) vs t with Ising crit- 
ical exponents /3 = |, v = 1, and z = 2.17 for A = 0.379 10, 
0.379 15, and 0.379 20 (from top to bottom) on a semiloga- 
rithmic scale. 



III. NUMERICAL ANALYSIS OF THE 
ORDER-DISORDER TRANSITION 

In this section, we present the simulation results, fo- 
cusing on the order-disorder transition. Rather than 
studying the Binder cumulant, we analyze how m{t) ap- 
proaches the steady state value. This approach is also 
known as the nonequilibrium relaxation method 1241. 

At criticality, the SM is expected to decay as [25, 26] 



m(t) - t-M"*), 
with the critical exponents /3, v, and z defined as, 



(12) 



91t~(A c -A)^ £~|A C -A|-", 



|A C -A|— , (13) 



where A c is the critical point, £ is the correlation length, 
and t is the relaxation time. For the two dimensional 
Ising model, (3 = g and v = 1 are known exactly (see, 
for instance, Ref. [27|), but z ~ 2.17 [26| is known only 
numerically; nonetheless it serves well for our purpose. 
In simulations, we set mo = 0.9 (initial condition) and 
observed how (M(t, L = 2 11 )) behaves up to t — 5 x 10 5 . 
When we study the finite-size scaling later, the finite size 
effect is argued to be negligible up to the observation 
time in this case, so {M(t,L = 2 11 )) can be regarded as 
m(t). 

In Fig. [21 we depict the behavior of m{t)t^^ vz ^ for 
three different values of A as a function of t on a semilog- 
arithmic scale, using the critical exponents of the two- 
dimensional Ising model. The numbers of independent 
simulation runs for A = 0.3791, 0.379 15, and 0.3792 are 
280, 912, and 160, respectively. In the ordered (disor- 
dered) phase, the curve is expected to veer up (down), 
and at criticality the curve should be flat, if the correct 
exponents are used. Thus, Fig. [2] supports the idea that 
the order-disorder transition in the 2DIM model is of the 
Ising type with the critical point A c = 0.379 15(5), where 
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FIG. 3. (Color online) Plots of (M(t, L))L Plv vs t/L z with 
Ising critical exponents at criticality on a double logarithmic 
scale. Inset: Plots of (\M(t, L)\)L p/v vs t/L z at A = A c . 
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FIG. 4. (Color online) Scaling collapse plot of m(t)h~ l ^ s as a 
function of th ' 16 on a double logarithmic scale with the Ising 
critical exponent 5 = 15. 



the number in parentheses indicates the error of the last 
digit. 

Note that the power-law behavior of m(t) is observable 
only for t > 10 4 , which implies that the 2DIM model has 
stronger corrections to scaling than the two-dimensional 
Ising model (for example, Fig. 1 in Ref. [26] shows that 
the two-dimensional Ising model is already in the scaling 
regime from t = 10). As we will see later, the strong cor- 
rections to scaling also plague the behavior of the Binder 
cumulant, which will be given as the reason why previous 
studies could not successfully report the universal value 
of the Binder cumulant (see Sec. HVBl) . 

To have further support for the Ising critical behav- 
ior, we also studied the finite size scaling. At criticality, 
scaling collapses for the magnetization M. and for the ab- 
solute value of the magnetization \M\ are expected with 
the scaling forms 



at criticality. Like the Ising model, the SM at the critical 
point is supposed to behave as [26 1 



(M(t,L)} = L-e/"f(t/L*), 
(\M(t,L)\)=L-P/»g(t/L% 



(14) 
(15) 



where / and g are (universal) scaling functions. To check 
the finite-size scaling at criticality, we simulated the sys- 
tems with sizes of L = 2 7 , 2 s , 2 9 , and 2 10 . The numbers 
of independent simulation runs for L — 2 7 , 2 s , 2 9 , and 
2 10 are 160 000, 40 000, 10 000, and 2504, respectively. 
The resulting scaling collapse is presented in Fig. [3] which 
indeed shows a nice collapse of (M(t, L)) onto a single 
curve when the Ising critical exponents are employed. 
The scaling collapse is also nice for the average of the 
absolute value of M. (inset of Fig. [3]). 

As a by-product of the finite-size scaling, we can esti- 
mate the time after which the finite-size effect becomes 
significant to be t « 0.1L Z (about 1.5 x 10 6 for L = 2 11 ) 
and, in turn, we affirm that the relaxation dynamics of 
m(t) is faithfully presented in Fig. [3J 

Furthermore, we also studied how the symmetry- 
breaking field h affects the behavior of the magnetization 



m{t) = h^ s H(th vz/ ^), 



(16) 



where 5 = 15 for the two-dimensional Ising model [271 ] 
and H is a universal scaling function. It is worth while 
to investigate a scaling collapse of m(t)/i -1 /' 5 plotted as 
a function of th uz ^^ s \ using the Ising critical exponents. 
Note that for the two-dimensional Ising model vz / ([38) ~ 
1.16. 

To introduce the symmetry-breaking field, we follow 
the idea in Ref. [28[. Now the transition rates take the 
form 

/ 1 - {-IY+3 \ 

W C „C =1 hr^h <$„ n*„ n 



+ S a lc ,1^2'^2-r S a v fiSn y ,k, (17) 



k 

y k=l 

where i and j are components of the lattice vector x 
and < h < 1. Recall that we have set fcAfc = A for 
k = 1,2, 3. By Eq. (fl"7|) . adsorption on the sublattice E is 
more probable than on the sublattice O, which eventually 
breaks the symmetry between even and odd absorbing 
states. 

With these modified transition rates, we simulated a 
system with size L = 2 9 at A = A c for different values 
of h. In actual simulations, we have only to change the 
probability of adsorption on a vacant site in the sublattice 
O without occupied nearest neighbors to (1 — h)At. The 
resulting scaling plot is depicted in Fig. [4] which shows 
a nice scaling collapse. It again supports the idea that 
the 2DIM model should belong to the Ising universality 
class. To make sure that the finite-size effect is negligible, 
we also studied systems with size L = 2 8 and obtained 
almost same figure as Fig. 2] (not shown here). 
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FIG. 5. (Color online) Plots of p(t) vs t~ 1/z with the Ising 
dynamic exponent z — 2.17 at criticality. The straight line is 
the result of a linear function fitting. Inset: Plots of (8p) 2 vs 
t at A = A c on a semilogarithmic scale. 



IV. DISCUSSION AND SUMMARY 
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FIG. 6. (Color online) Plots of U(t,L) vs i at criticality 
(A = A c ) on a semilogarithmic scale (L = 2 7 , 2 8 , 2 9 , and 
2 10 from left to right). As a guide to the eyes, we also plot a 
straight line indicating the universal Binder cumulant of the 
Ising universality class, 0.611. Inset: Plots of U(t,L) vs t in 
the disordered phase (at A = 0.3795 > A c ) for L = 2 6 , 2 7 , 2 s , 
and 2 9 from left to right. 



Up to now, we have shown that the critical behavior of 
the order-disorder transition in the 2DIM model is of the 
Ising class. In this section, we will discuss the behavior of 
the density of active sites p(t) and the Binder cumulant 
at the order-disorder transition point A c , which will be 
compared with the similar studies in Ref. [l9[. Also, to 
confirm the universality, we discuss the critical behavior 
of the absorbing phase transition. 



A. Behavior of p and its fluctuation at A = A c 

In Ref. [l9|, diverging fluctuation of the order parame- 
ter of the absorbing phase transition was suspected to be 
a possible reason why the two dimensional IMD model 
should not belong to the Ising class. However, this or- 
der parameter, in our case <j> defined in Eq. (Q}, seems 
related to the energy density of the Ising model in that 4> 
is measured as a correlation between nearest neighbors. 
Since the fluctuation of the (Ising) energy is the specific 
heat which diverges logarithmically at criticality in two 
dimensions, it is actually plausible that the fluctuation of 
4> (times system size) diverges at criticality even though 
the 2DIM model belongs to the Ising class. 

To confirm that pit) = (4>) indeed is linked to the (av- 
erage) energy density of the Ising model, we analyze how 
pit) behaves at the order-disorder transition point. Since 
the energy at criticality approaches the steady-state value 
in a power-law fashion with exponent [vd — l)/(i/z) [2fj| . 
pit) is expected, if it is indeed related to the energy, to 
approach the steady-state value p* in such a way that 



-pit) »t-l/» 



(18) 



where we have set v = 1 (the Ising critical exponent) 



and d = 2. Hence if pit) is plotted as a function of 
t- x l\ the curve becomes straight for small t x l z and 
approaches p* as t~ x / z — >• (equivalently, t —> oo). As 
Fig. [5] reveals, p approaches the ordinate as a straight 
line for i" 1 ' 2 < 0.02, as anticipated. Note that the time 
i -1 / 2 = 0.02 roughly corresponds to t — 5 x 10 3 after 
which m(t) enters the scaling regime (see Fig. [2]). 

We also analyzed how the fluctuation of the active site 
density defined as 



i5pf= Urn L 2 ((0(t,i) 2 )-(0(t,i)) 2 



(19) 



behaves at A = A c . The inset of Fig. [5] shows logarithmic 
behavior of (<5p) 2 as in the two-dimensional Ising model. 

Thus, we conclude that the active site density 4>(t, L) 
is indeed associated with the energy of the Ising model. 
Actually, the logarithmic behavior of (<5p) 2 is compatible 
with the slow divergence of the fluctuation observed in 
Ref. [H. 



B. Binder cumulant 

Since the Binder cumulant is believed to take a uni- 
versal number at criticality, we should study whether the 
Binder cumulant at A = A c approaches to the universal 
value as L — > oo. Denning the Binder cumulant at time 
t as 



U(t,L) = 1 - 



3(Mit,L)*y 



(20) 



we numerically study how Uit, L) behaves for different 
values of L. 
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FIG. 7. (Color online) Plots of p(t)t 0A5 vs t near the critical 
point of the absorbing phase transition on a semilogarithmic 
scale. 



In Fig. [6] we present simulation results for U(t, L) for 
A = A c = 0.379 15 and for A = 0.3795 > A c (Inset). 
For A = 0.3795, the numbers of independent samples 
simulated for L = 2 6 , 2 7 , 2 8 , and 2 9 are 200 000, 50000, 
14000, and 4000, respectively and data for A = A c were 
collected while we studied the finite size scaling. 

Since the system has absorbing states and any finite 
system will eventually fall into one of the absorbing states 
even in the active phase, there are obviously two char- 
acteristic time scales. One is r q when the system enters 
the quasi-stationary state and the other is r Q when the 
system falls into one of the absorbing states. At A = A c , 
T q diverges with system size as r q ~ L z , but r Q should 
increase exponentially with L because the SB transition 
point is in the active phase of the absorbing phase tran- 
sition. Hence, to find the universal value of the Binder 
cumulant at the SB transition point, the observation time 
should be larger than r q but much smaller than r a . Actu- 
ally, except for the case of L = 2 6 , no simulation results 
in an absorbing state and, even for L = 2 6 , only < 0.2% 
of simulation runs falls into an absorbing state up to the 
observation time. Hence, in our analysis, the Binder cu- 
mulant is not influenced by the existence of absorbing 
states. 

At A = A c , U(t, L) in the (quasi-)stationary state in- 
crease with system size but shows a clear signature of 
saturating behavior to the universal number 0.611 (29| 
as I — > oo. Note that if the system size is not large 
enough, the Binder cumulant could miidetify the critical 
point. The unexpected behavior of the Binder cumulant 
should be attributed to the strong corrections to scaling 
already observed in Fig. [2] 

The inset of Fig. |6] depicts the behavior of U (t, L) at 
A = 0.3795 > A c (disordered phase). If the system size is 
not larger than L = 2 s , one may conclude that the criti- 
cal point is around 0.3795, with the value of the Binder 
cumulant around 0.59, which is comparable to the value 
reported in Ref. [19j . Hence we conclude that the critical 



point reported in Ref. [l9j is actually in the disordered 
phase. 



C. Absorbing phase transition 

Finally, we discuss the critical behavior of the absorb- 
ing phase transition. Since the symmetry is already bro- 
ken, it is expected that the model should belong to the 
DP class [la, Hit- To confirm this, we study the system 
with the initial SM m = 0.1. In Fig. [3 we plot p(t)t 0A5 
as a function of t on a semilogarithmic scale, where 0.45 is 
the critical exponent of the DP class. For A = 0.366 750, 
the curve becomes flat from around t = 5 x 10 4 . In the 
active (absorbing) phase, the curves veer up (down) as 
usual. Thus we conclude that the critical point of the 
absorbing transition is A Q = 0.366 750(5) and the critical 
behavior is of the DP class. 

Note that exponential decay of p in the absorbing 
phase is observed in Fig. [3 which might look inconsis- 
tent with the power-law decay in the whole absorbing 
phase reported in Ref. 19]. However, there is a clear 
distinction. Since the initial SM is nonzero in our case, 
coarsening has not played any role. Indeed, we also ob- 
serve power-law behavior in the absorbing phase if tuq is 
set to zero just as in Ref. [l9| (data not shown). 



D. Summary 

To sum up, we studied a model of two-dimensional in- 
teracting monomers, focusing on the order-disorder phase 
transition. Numerical analysis showed that the 2DIM 
model should belong to the Ising universality class, con- 
trary to a recent claim We observed that analysis 
of the Binder cumulant is not an efficient method to find 
the critical point in two dimensional models with two 
symmetric absorbing states. We also reconfirmed that 
the absorbing phase transition occurring in the 2DIM 
model after the symmetry is broken is described by two- 
dimensional directed percolation. 

Although we did not directly study the interacting 
monomer-dimer model, we believe that the conclusion in 
this paper should be applicable to the IMD model studied 
in Ref. [I9j because of the universality hypothesis. Since 
the two different models, IMD and 2DIM, have strong 
corrections to scaling at the symmetry-breaking transi- 
tion point, unlike the Ising model, the origin of these 
strong corrections seems to be related to the presence of 
an absorbing state even for A > A c (disordered phase). 
If this is the case, it is an interesting question as to why 
and how the absorbing states affect the corrections to 
scaling; this is beyond the scope of the present paper 
and is deferred to a later presentation. 
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